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Abstract 

Using AdS/CFT, we compute the Fourier space profile of (trF 2 ) generated by a heavy 
quark moving through a thermal plasma of strongly coupled M = 4 super- Yang- Mills theory. 
We find evidence of directional emission from the quark whose description includes gauge 
fields with large momenta. We comment on the possible relevance of our results to relativistic 
heavy ion collisions. 
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1 Introduction 



In PUI2], a classical solution of string theory is described that is dual in the sense of AdS/CFT 
[HJ IU E] to an external quark passing through a thermal plasma of M = 4 super- Yang-Mills 
theory at large N and strong 't Hooft coupling gy M N . The string is treated in the test string 
approximation: its back-reaction on the geometry is not considered. The string dangles into 
AdS^-Schwarzschild from an external quark on the boundary which is constrained to move 
with constant velocity. The string trails out behind the quark and exerts a drag force 

<^P = ^V9ym n T 2 v _ 

dt 2 % /T^2 

on the quark. Here v is the speed of the quark, and T is the temperature of the plasma, or 
equivalently the Hawking temperature of the horizon of AdSs-Schwarzschild. The diffusion 
constant D = 2/ (nT ^g YM N) implied by ((TJ) was derived independently in 0, also using 
AdS/CFT. 

In the gauge theory, energy loss results from gluons (or superpartners of gluons) radiating 
off the heavy quark and interacting with the plasma. We should ask: How energetic are these 
radiated gluons? At what angle do they come off relative to the velocity of the heavy quark? 
To the extent that such questions can be posed in a gauge-invariant manner, AdS/CFT 
should be able to provide an answer. The aim of the present paper is to shed some light 
on these questions by computing the profile of (trF 2 ) in the boundary gauge theory. To do 
this we compute the linear response of the dilaton field to the string, which is a first step in 
computing its back-reaction on the /WS' 5 -Schwarzschild background. Actually, what we will 
extract in the end is the vacuum expectation value (VEV) of the operator in M = 4 super- 
Yang-Mills which couples to the dilaton. This is not quite trF 2 , but rather the lagrangian 
density plus a total derivative: in mostly plus signature, 

F 2 = tr (-F^ n + 2X I D 2 m X I - 2i\ a a m D m \ a + more interactions) , (2) 

where X 1 are the six adjoint scalars, A a are the four Weyl adjoint fermions, D m is the 
gauge-covariant derivative, and a m = (—1, — (?) where a are the Pauli matrices. 

The near field of the heavy quark is just the Coulomb color-electric flux, appropriately 
Lorentz boosted. The contribution of this near field to (Op2) can be computed analytically, 
following [7], and it has nothing to do with energy loss. When it is subtracted away from 
(Ops), the remainder is peaked at momenta many times larger than the temperature. The 
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information in (Opt) is complementary to in that it helps identify the energy scale 
at which dissipative phenomena occur but does not so clearly indicate the overall rate of 
dissipation. More complete information could be extracted from (T^), which could also 
be computed via AdS/CFT but requires a more technically involved treatment of metric 
perturbations. 

Several related papers (HI El E3 IH] appeared recently, all aiming to describe at some 
level energy dissipation from a fundamental quark into a thermal plasma using AdS/CFT. 
The interest in this topic owes to a possible connection with relativistic heavy ion physics. A 
distinctive feature of RHIC experiments f21EEElE3 is jet-quenching, which is understood 
as strong energy loss as a high-energy parton passes through the quark-gluon plasma formed 
in a gold-on-gold collision. 

The organization of the rest of this paper is as follows. In section|2]we explain the classical 
supergravity calculation that leads to (Op?)- Similar calculations were carried out in ^H] 
for a string undergoing small oscillations around certain static configurations in AdS§. All 
the supergravity computations are done in five dimensions, but the final answer is the gauge 
theory quantity (Of 2 ) as a function of the coordinates (t, x 1 ,^ 2 ,^ 3 ) of Minkowski space. 
(Actually we will find it easier to pass to momentum space early in the computation.) One 
must solve a boundary value problem in order to extract (Ops). Numerical techniques for 
doing so and results for several different choices of v are described in sectional We conclude 
in section 0] with a discussion of the possible relevance of our work to recent experimental 
results. 

2 Dilaton perturbations 

The background geometry is the well-known ArfS's-Schwarzschild solution, 




(3) 



and useful relations include 
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In static gauge, the string worldsheet is described as 



X»(t,z) = (t X\t,z) z) 
X 1 {t,z)=vt + £(z) £{z) = 



z H v ( 1-iz/zh , l + z/z H 

log ; — h % log 



(5) 



Ai \ ° 1 + iz/zn 1 — z/zh , 

To compute the dilaton response to this string, one starts with the following action: 

1 



S = J d 5 x V^G 



4«1 



(dcpy 



[ d 2 ae^ 2 V=9 g aP = G^d a X' i dpX v . (6) 

27TQ!' J M 



Here a and (5 refer to worldsheet coordinates a a = (r, a), and 
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(7) 



where G5 is the five-dimensional gravitational constant. To derive the dilaton equation of 
motion, it helps first to rewrite the whole action as a single volume integral (we refrain briefly 
from choosing static gauge): 



S = d 5 x V^G 
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-2 1 / d 2 ae <P/2 V_9_ A 5 



{d<j>y - 



2ira' 



5 b (x» - X»(cr)) 



The five-dimensional delta function in (jHJ) is a product of standard Dirac delta functions. 
So, for instance, 

5 5 (x^) = S(t)S(x l )S(x 2 )S(x 3 )S(z) . (9) 



The linearized equation of motion can now be straightforwardly derived as 

=<S 5 (x" - X"(a)) . (10) 



□ 



duV^GG^drf = 3=-^-. I d 2 a V 9 



V^G 



2na' 



V^G 



and by passing to static gauge one may explicitly perform then the remaining integral in 



j 



K 5 V 9 r/ 1 V l 



Six 1 -X\t,z))5(x 2 )5(x s ). 



11] 
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In the spirit of finding the steady-state, late-time behavior, we assume that <fi depends 
on x 1 and t only through the combination x 1 — vt. After computing 



V-9 = z_ 
V^G L 3 



vT 



(12) 
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one can easily show that □ = J simplifies to 



z 3 d z -d z 



-^)d! + d* + dt 



V 2 z 



2vra' L 



Six 1 -vt-£{z))5(x 2 )5(x 3 ). (13) 



This partial differential equation can be attacked by Fourier transforming: 

d 3 k 



(j)(t,x, z) 



and similarly for J. Then one has 



D iki (x 1 — vt)+ik2X 2 +ik3X 3 



(14) 
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where k\ = k\ + k\. All dimensionful factors drop out of the differential equation when we 
introduce rescaled variables 



Ki = z H ki K± = z H k± y 



Zr 



2na'L 1 
j k{v) = -ZT3 n o Mz) ■ (16) 
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Then h = 1 — y A and 



h 



fd y - dy -{l--)Ki-Kj 



1 + iy 



vKt/4 



1) 



ivKi /4 



i-yj 



(17) 

There doesn't appear to be a solution to (fT7|) in terms of known special functions. However 
it can be solved in two interesting limiting regimes: 

• Near the horizon, y is slightly less than 1, a better choice of radial variable is 
Y = log(l — y). The leading terms in the differential equation near the horizon (that 
is, for large negative Y) are 



81 + ' vKl 
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± e Y e -ivK 1 (Y+w/2-log2)/4 



which is also the equation of motion for a simple harmonic oscillator with a complex 
driving force. The solutions are 



e y /4 
1 - WKJ2 1 



-ivK 1 (Y+7r/2-lo g 2)/4 _|_ Q+ e ivKiY/4 + 
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where C K are arbitrary constants. The standard boundary condition at a black hole 
horizon is to choose a purely infalling solution. This means that in the near-horizon 
limit, <p should depend on t and Y only through the combination t + ZhY/4, not 
t — zhY/A: the quantity ZhY/4 is essentially the tortoise coordinate. Thus C K = 0. 

Near the boundary of AdS's-Schwarzschild, the leading terms in the differential 
equation are 

y 3 d y —d y 4> K = y , (20) 



and the solutions are 



W = -y +Ak + B K y\ (21) 



where Ak and Bk are arbitrary constants. Ak should be set to zero because there is 
no deformation of the lagrangian. Bk is proportional to (Op2). 

It is worth noting that the relation of Bk to (Of 2 ) involves a subtraction of contact 
terms. Conventionally, it is understood that 



(0 F2 (t,x)) = -i^liml^(t,f, z), 

Zftg 2^0 Z" 



but in the present case, the limit doesn't exist because of the y 3 term in ()21|). Fortunately, 
this term has no K dependence. Thus when passing back to real space, it is proportional to 
a delta function supported at the location of the quark. This delta function has an infinite 
coefficient, but if it is subtracted, the remaining contribution to (Op2(t,x)) indeed comes 
from Bk, and it is finite. The subtraction prescription has some arbitrariness: one could 
subtract off any finite multiple of the delta function at the same time, which corresponds to 
subtracting a i^-independent quantity from every Bk- 
Combining (gj), flUD, (EH), and fl22J), one finds 

(O p2 (t,x)) = -ir 3 T^g$ M NVT^ J 0_ e i^(* 1 -vt)+^ 2 +^ 3 }/zH BK ^ ( 23) 

In sectional we will quote results in units where zh = 1: this corresponds to T — l/V. 

For a wide range of K\ and the dominant contribution to Bk comes from the near 
field of the quark, which in position space is proportional to l/|x| 4 in the rest frame of the 
quark. Consider first the case v = 0. Following consider a string dangling straight down 
in AdS§. One obtains 
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This calculation is done in the absence of a horizon, or equivalently at zero temperature. 
Fourier transforming (J24)) leads to 



We have expressed the result in terms of the dimensionless variables (jlfij) with zh = 1/irT 
finite, even though T = physically. This is a bookkeeping trick to obtain a form that can 
easily be compared with AdSs-Schwarzschild results. 

For v 7^ 0, one may apply a Lorentz boost to the AdS§ string configuration considered 
in the previous paragraph. This describes an external quark moving through the vacuum at 
speed v. The result for _B]*? ar field in this case is 



This is the analytic form that we will subtract from numerically evaluated Bx to excise the 
near field but leave behind all the dissipative dynamics. 

3 Numerical algorithms and results 

The boundary value problem described in and below (|17j) is reminiscent of both the glueball 
calculations initiated in arid of quasi-normal modes in AdSs-Schwarzschild ^Hj- But 

there is an additional simplifying feature: all the equations are affine in <px — that is, they are 
linear combinations of 4>k(h), its derivatives, and functions of y that do not involve 4>k{v)- 
To see this, consider the following formulation of the horizon boundary condition. One first 
expresses the asymptotic solutions (fi n c&r,K and 0f ar ,i^ as a sum of the inhomogenous solution 
and the permitted homogenous solution. Explicitly, for the near-horizon solution, 




(25) 




(26) 



near, K 




ivK 1 (Y+n/2-log2)/4: 



(27) 



'near, P, K 



1 - ivK x j2 e 



The Wronskian 



Wnear(2/) = (<j> K (y) ~ 0ncar,P,^(l/))0 ne ar,H,^(l/) ~ WM ~ 0near,P,if iv) )0near,H,i^ (y) (28) 
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is a measure of how close the numerically computed function 4>k{u) is to the analytic ap- 
proximation <fi near ,K- Because the horizon is a singular point of the differential equation, one 
must impose the boundary condition W neai (yi) = at a point y± slightly less than 1, which 
is to say slightly outside the horizon. The quantity W n ear(z/i) is indeed a linear combina- 
tion of 4>K(y), $K{y)i an d a (^-independent function known in terms of <^ n ear,H,A:(2/) an d 
<t>mmxp,K{y)- One ma y similarly formulate a boundary condition Wfax(yo) — which is also 
affine in <f>K- The point yo should be chosen slightly greater than 0, which is to say close to 
the boundary of Ac/Ss-Schwarzschild. 

There are special methods to solve boundary value problems of the type just described, 
where both the differential equation and the boundary conditions are affine, which are more 
efficient than standard shooting algorithms. Mathematical NDSolve incorporates such meth- 
ods internally [201 • But we have found that we achieve greater numerical accuracy using a 
home-grown shooting method where Bk is guessed and then adjusted to make C K = 0. 
Accuracy was further improved by finding power series corrections to the asymptotic forms 
()19|) and (J21J) . A satisfactory choice of cutoff points was yo = 0.01 and y\ = 0.99. The 
numerical challenge increases as K\ and K± increase, requiring more CPU time. As we will 
see in figure El B% is significantly weighted toward K larger than 10 when an appropriate 
phase space factor is included. So it would be worthwhile to have some alternative method 
adapted to this regime, perhaps based on a WKB approximation. 

We take advantage of the axial symmetry of the problem to express Bk = B(Ki, K±) 
where K± = y 'K\ + K\. Because (j)(t,x,z) and (Of^,x)) are real, it must be that 
B(—Ki,Kj_) = B(Ki,K±)*. It is easy to see that this condition is enforced by the dif- 
ferential equation. Our results for B(Ki, K±), with the near field ()26|) subtracted, are shown 
in figures ^ and 121 A good match to the near field form f)26|) was obtained: for K± > 10 
the deviations are at the level of tenths of a percent. These deviations are interesting and 
can be seen in magnified form in panes b, d, f, and h of figure 121 Much of our discussion in 
section 0] will hinge on these high- momentum tails. 

4 Discussion 

Before attempting a comparison of our results with recent RHIC results, we will give a brief 
summary of how the measurements of interest are done. The reader is warned that we are 
non-experts and is referred to the experimental literature — for example [2H 1221 12H] — for an 
authorative account. 
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a) Re B(Ki,K p ) for v=0.75 



b) -Im B(Ki,K p ) for v=0.75 




Figure 1: Contour plots of the real part and minus the imaginary part of B(Ki, Kj_) for 
several values of v. The near field contribution (|26|) has been subtracted. B(K\, K±) is 
proportional to the K-th Fourier mode of {O f %)\ see f)23|) . In each plot, the white region is 
closest to zero, and the black region is the most positive. 
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a) |B(Ki,K p ) | for v=0.75 b) K p |B(Ki,K p ) | for v=0.75 




10 20 30 40 10 20 30 40 



Kj Ki 

Figure 2: The absolute value of B(Ki, K±) with and without the phase space factor K±. 
The near field contribution ()26|) has been subtracted. The green dot is the recoil energy of a 
thermal gluon: see (|3UJ1 . The dashed red lines indicate the direction in which K±\B(Ki, K±)\ 
is largest: see the discussion around In each plot, the white region is closest to zero, 

and the black region is the most positive. 
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Consider the following scenario: 

1. Two highly energetic partons collide near the surface of the hot dense matter produced 
in a relativistic heavy ion collision. After the collision, the partons have large transverse 
momentum. 

2. One parton escapes without interacting significantly with the quark-gluon plasma 
(QGP) and fragments in vacuum into what is termed the near side jet. 

3. The other parton travels through the QGP. Its evolution into observed particles is 
strongly affected by its interaction with the QGP. If it weren't for these interactions, 
this parton would simply fragment into an away side jet, approximately back-to-back 
with the near side jet. 

Because of difficulties in unambiguously identifying jets, a standard strategy is to look for 
angular correlations between two energetic charged particles: the trigger particle, which is 
presumed to be part of the near-side jet, and the partner particle, which is the putative 
probe of jet-quenching. Histograms of the azimuthal angle A0 between these two particles 
invariably show a peak at small angles, which means that the partner particle is often 
part of the near-side jet. A peak at A0 = n is evidence for an away side jet. In central 
collisions, the peak at A0 = ir disappears [2J or even splits [221 122j In [22], the trigger 
particle is required to have 2.5GeV/c < < 4.0GeV/c while the partner particles has 
l.OGeV/c < pt < 2.5GeV/c, and for central collisions a broad peak is observed roughly 
between A</> = 1.6 and A0 = 2.6. (All angles will be quoted in radians.) There is actually a 
minimum at A<p = tt. 

The recent theoretical literature on jet-quenching, with which we have less familiarity 
than we would like, offers several possibilities. Among them are scenarios [2U [2Sj where 
the QGP affects fragmentation by recombination of thermal quarks with the parton shower; 
extensions of traditional QCD methods such as the twist expansion j2H]; predictions of a 
coherent high momentum ridge of color flux emanating from the quark fZT\ |2*H] ; and related 
discussions of a QCD "sonic boom" giving rise to conical collective flow [2*^113*0] . 

In the backdrop of these experimental and theoretical investigations, it is interesting to 
say what we can about the energy flow and spectrum of particles radiated from the heavy 
quark described in the previous sections. The hazards of comparing strongly coupled M = 4 
super- Yang-Mills with real-world QCD are well known: for a brief summary, see j2J- To these 
difficulties we must add that we have treated the quark as infinitely massive, whereas the 
experimental results we have referred to do not include heavy-quark tagging. Also, it would 
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be better to know (T^ u ) in addition to (0_f 2 ) : energy flow is most crisply captured in the 
Poynting vector Si = T Qi . Finally, it would be desirable to go to larger K\ and K±, which 
requires either CPU-intensive numerics or an improved calculational method, as discussed 
near the end of section El 

Objects deep inside AdS$ are understood to correspond to soft field configurations in the 
dual CFT, while objects near the boundary correspond to more localized configurations. So 
a reasonable expectation based on ^ |2] is that the profile of (Of 2 ) would have the form 
of a wake, consistent with the ideas of [23 I2H I2H ED]- The Fourier space profiles shown 
in figures H] and |21 suggest a slightly different, possibly complementary picture. It helps our 
intuition to use explicit numbers. Let's set 

T = — GeV = 318 MeV . (29) 

7T 

This is in the upper range of temperatures for the QGP, and it is a convenient choice for 
us because the K\ and K± axes in figure ^ and 121 can then be read in units of GeV/c. 
Another interesting number is the typical final energy of a free massless particle that collides 
elastically with the heavy quark. To compute this we take the initial momentum of the 
massless particle to be of magnitude T and directed perpendicular to the heavy quark's 
velocity. If the perpendicular component of the massless particle's momentum doesn't change 
during the collision, then its final energy is 

1 + v 2 

E f = -^^T = 6.2 GeV for v = 0.95. (30) 
1 - v 2 

We have indicated Ef for the various velocities with the green dots in panes a, c, e, and g of 
figure 121 If the gauge theory were almost free instead of strongly coupled, we would expect 
the energy loss to be dominated by collisions of the type that led to (J30|) . 

For v = 0.95, \B(Ki, K±)\ is peaked in a range of momenta between 2 and 7GeV/c (the 
black region in figure |2J). Because Opi ~ trF 2 starts with bilinears in the fundamental 
fields, this would correspond to radiated particles with momenta between 1 and 3.5GeV/c: 
less than the Ef of (J30|) by a factor of a few. For v = 0.99, half the momentum at which 
\B(Ki, K±) \ is peaked is less than Ef by a similar factor. These considerations encourage the 
view that dissipative events involve several quanta interacting with each other as they recoil 
from the heavy quark. This is broadly consistent with the picture of a coherent co-moving 
high momentum ridge dissipating energy from the heavy quark. But as we will see below, 
the peak regions of \B(K\, K±)\ may not dominate the dissipative physics. 
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Panes b, d, f, and h of figure El show that if one multiplies \B(K\, K±) \ by the factor K± 
that would arise in an integration over momentum space, the result is directionally peaked. 
This again brings to mind the picture of dissipation through radiation carried mostly in the 
high momentum ridge. The opening angle 9 between the heavy quark's velocity and the 
directional peak of K±\B(K 1 , K±)\ depends strongly on the speed: 
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0.41 
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0.17 



The values of 9 in ()31|) were determined by setting the K± derivative of K±\B(Ki, K±)\ to 
zero at fixed and large Ki, then taking the appropriate arctangent function to find 9. 

The phase space factor K± makes an enormous difference to the dominant momentum 
scale. In K±\B(K\, K±)\, momenta many times Ej dominate. Indeed, along the preferred 
direction, K±\B(Ki, K±) | seems to level off at a finite value as K increases. Evidently we have 
not explored sufficiently high momenta to discern whether the region where K±\B(K\, K±)\ 
is above a finite threshhold has finite volume. 1 

To recap: the plots of K±\B(K\, K±)\ not only indicate directionality, but also suggest 
that highly energetic fields are an important part of the description of the radiation process. 
To appreciate just how energetic, note that a charm quark moving in vacuum with v = 0.95 
has energy 4.5 GeV, while a b quark with this speed has energy 15 GeV. If K±\B(Ki, K±)\ 
can be used as an approximate guide to the spectrum of radiated particles, the single particle 
energy could easily be in the 10 GeV ballpark. Recoil would obviously become an important 
consideration if a real-world c or b quark emitted a particle even approaching this range. 
This would substantially increase the opening angle 9. And it would encourage the idea that 
the QGP enhances fragmentation processes at energies close to the kinematic limit. 

There are two main reasons to treat with particular caution a "prediction" from AdS / CFT 
that heavy quarks should undergo fragmentation near the kinematic limit: 

1. We have not made a quantitatively precise connection between (Ops) in Fourier space 

and the spectrum of radiated particles. Indeed, the peak region of B(Ki, K±) and 

its high-momentum tails send conflicting messages about the spectrum. We believe 

the tails are important, but it may be that they have to do mostly with fields near 

1 Note that in the large momentum region of the plots shown, we are subtracting a quantity, £?™ ar fiold , 
which scales linearly with momenta. The remainder, B(K\,K±), scales roughly as 1/K in the region in 
question. This evidently requires substantial numerical precision. All internal checks of our numerical 
results suggest they are robust, but the importance of large K tails to our discussion is the reason we say it 
would be value to have WKB methods in hand. 
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the quark rather than radiative dynamics. The question of the spectrum of radiated 
particles should be revisited purely within the context of AdS / CFT with the VEV of 
the stress-energy tensor in hand, and preferably with semi-analytic methods to buttress 
numerical analysis of the high- momentum tails. 2 

2. Relating hard processes in strongly coupled M = 4 super- Yang-Mills and QCD is es- 
pecially perilous. Elementary scattering processes with large momentum transfer can 
be treated perturbatively in QCD. In strongly coupled M = 4 super- Yang- Mills the 
general expectation is that they cannot. But one should bear in mind that many am- 
plitudes of M = 4 are protected against all loop corrections. It would be interesting 
to inquire whether amplitudes for gluons scattering off an external quark have non- 
renormalization properties. This discussion recalls the basic conundrum of the connec- 
tion between AdS/CFT and RHIC: are near-extremal D3-branes merely an analogous 
system to the QGP, or can they capture the dynamics of real-world QCD above the 
confinement transition sufficiently precisely to be a useful guide to RHIC physics? 

Fragmentation near the kinematic limit seems to us consistent with the broad peak in 
A(f> observed in [22] • But the energy ranges for the hadrons in are substantially lower, 
relative to the temperature, than the energies we have discussed in relation to AdS/CFT. 
Recall that the upper limit on of the partner particle is 2.5 GeV/c. If the typical energy 
of the partner particles is sufficiently low, it would be a blow to the picture of enhanced 
high-energy fragmentation. Of course, without tagging most of the partons studied in [22] 
may be presumed to be light quarks or gluons. 

In summary, the calculations we perform are based on the trailing string picture of P}|2], 
which naively supports the notion of a coherent wake of color fields with the heavy quark at 
its tip. We do find evidence for a directional "prow," which becomes more and more forward 
as the speed increases. It seems that a full description of this prow involves high-momentum 
gauge fields. This may be a hint that, with a realistic cutoff on the quark mass imposed by 
hand, the quark could be deflected significantly by a single radiative event. 

The drag force computed in 13 is a time-averaged quantity which provides no 

direct information about the energy scale of radiated particles. Calculating color-singlet 

VEV's in the boundary theory gives considerably more detailed information. Despite the 

hurdles string theory faces in connecting to relativistic heavy ion collisions, we hope that the 

trailing string picture can be further exploited to understand energy loss in the QGP. 

2 Indeed, a computation of the stress tensor gives clear-cut evidence at smaller wave-numbers for a wake 
in sense usually meant by phenomenologists |31| . 
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